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ABSTRACT 


Different agencies using the tracking data from the GEOS geodetic tracking 
systems employ different mathematical models of the atmosphere or different approx- 
imations to the same model. Consequently, different refraction corrections are ap- 
plied by different users of the same data. This paper compares the magnitude, at 
different elevation angles, of some of the afferent tropospheric refraction correction 
formulations utilized by Goddard Space Flight Center, Wallops Island, Smithsonian 
Astrophysical Observatory, Army Map Service, Naval Air Systems Command, Central 
Radio Propagation Laboratory, and Eastern Test Range for refraction corrections to 
the laser, GRARR, SECOR, TRANET, and C-band radar measurements of range, 
range rate, and elevation angle of earth satellites. At a 15-degree elevation angle 

Q 

and for a nominal surface refractivity of N g = 313 x 10 , the maximum differences 

between the above formulations are about 2.3 meters in range, 1.4 cm/ sec in range 
rate, and 31 seconds of arc in elevation angle. 
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GEOS SATELLITE TRACKING CORRECTIONS FOR 
REFRACTION IN THL TROPOSPHERE 


INTRODUCTION 

The refraction correction formulas for correcting the elevation angle (E) , 
o 

range (R) , and range rate (R) tracking data from the GEOS geodetic tracking systems 
are not standardized. Different refraction correction formulas are applied by differ- 
ent users of the same data. The different formulas arise from the use of different 
mathematical models of the atmosphere or from different approximations to the same 
model. 

This document presents some of the different sets of tropospheric refrac- 
tion correction formulas now being used on GEOS data, and compares the magnitude of 

o 

the corrections to E, R, and R among the different formulations at different elevation 
angles. 

o 

Five sets of refraction correction equations for E, R, and R are being 

used by the Goddard Space Flight Center (GSFC) on GEOS data. These are found in 

a refraction study done for GSFC in 1964 by Dr. J. J. Freeman and in four GSFC 

orbit determination programs (ODP's). The four ODP's are the operational orbit 

differential correction (DC) program, the NONAME definitive orbit program, the 

geodetic data analysis program (GDAP), and the network analysis program (NAP). 

o 

These refraction formulas for E, R, and R are compared with each other in this 
document. The formulas for refraction correction in E and R are also compared 
with the Army Map Service (AMS) and Engineer Topographic Laboratory (ETL)SECOR 
R refraction correction formulas, with the Wallops Island C-band radar E and R re- 
fraction correction formulas, and with the GSFC and Smithsonian Astrophysical Ob- 
servatory (SAO) optical laser R refraction correction formulas. The formulas for 
o 

R are compared with the Applied Physics Laboratory (APL) Hopfield and Naval 
Weapons Laboratory (NWL) Astro range rate refraction correction formulas. 

All the E and R refraction correction formulas are also compared with 
the Central Radio Propagation Laboratory (CRPL) tables of the reference exponential 
troposphere (Reference I) and with the results of ray tracings through the same CRPL 
reference exponential troposphere using the Eastern Test Range (ETR) REEK program. 

In addition, the R refraction correction formulas are compared with values of = 

o dt 

E , where is the rate of change of REEK ray-trace range error (AR) 

with respect to E obtained by differentiating the polynomial fitted to the curve of 
o 

AR vs E, and E is the time rate of change E. 
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n 

< 4 > 

In Appendix A, the tropospheric refraction correction equations are first 
given as they appear in the referenced source documents and then are converted to a 
common notation and functional form to simplify the comparisons of these equations. 

o 

TROPOSPHERIC REFRACTION CORRECTION FORMULAS FOR E. R, AND R 


o 

The E, R, and R tropospheric refraction correction formulas used with 

GEOS data are listed in Tables 1,2, and 3, respectively. All formulas listed in these 

o 

tables employ nominal or first-order correction functions for E, R, and R, defined 
as: 


AE = N ctnE 
to s 

AR = HN cscE 
to s 

a/L = HN E ctnE cscE 
to s 


where N, = refractivity at station (nominally N = 0.000313) 
s s 

H = scale height of exponential model of troposphere from the 

expression N (h) = N e”^' H (nominally, if N = 0.000313 
s s 

at h = 0, then from the CRPL model, H = 6951.25 meters) 
E = measured elevation angle 

o o 

E = time rate of change of E. In these comparisons E is cal- 
culated from an overhead pass for a satellite at a height of 
(4/3) x 10 meters (see Appendix B) 


% 


t 



At low-elevation angles the nominal correction functions all approach in- 

o 

finity, whereas the true E, R, and R tropospheric radio refraction corrections at the 
horizon (E = 0°) are all finite with values of about 0.7°, 100 meters, and 200 cm/sec, 
respectively. Therefore, some of the formulations employ restraining functions 
which, when multiplied with the nominal correction functions, produce new correction 
functions to yield finite refraction corrections closer to the true values at low- 
elevation angles. 

o 

The limiting values for the E, R, and R tropospheric refraction correc- 
tions obtained from the formulas in Tables 1, 2, and 3 and from the CRPL tables are 
given in Table 4 for the horizon case (E = 0°) and for the zenith case (E = 90°). For 

the zenith ease, the E tropospheric radio refraction correction vanishes. At the 

o 

satellite point of closest approach (PCA) to the station, the R tropospheric radio 
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Refraction Appendix A Range Refraction Correction AR 

formulation Equation A R (meters) = Obs - Corr. 



fi 


where A R HN cscE IlAK. sccE r first order correction 



Table 3 

Range Rate Refraction Correction Equation 


Refraction 

Formulation 


Appendix A 
Equation 


Range Rate Refraction Correction, A R 
o 

£. R (metera/secl = Oba-Corr. 


GSFC DC 

A3a (E < 10°) 


GSFC DC 

A3b (E > 10°) 


GSFC Freeman 

A6 

A 8 to[> 

GSFC NONAME 

A9 

A»io[ 

GSFC GDAP 

A12a 

>K[ 


GSFC NAP-1 


A PL TRANET 


NWL TRANET 


n (1 ♦ 0.026 escE) J 

7200 2 ' 

R / o o 

J 1 + 0.0045154 cac^E + (1 + 0.0045154 csc*E)- 


.» [2.7432 

*° HN . / 

L VI * «• 

[««] 


2 2 
004 esc E + (1 ♦ 0.004 esc I 


AR ^ 
AR to HN S 


In equation A15, 


/ 2R \ r ~ 5 5 o - ^ + 

T'E) = 1 * -r slnE ] VH. ♦ 2R H + R ain^E - R sinE + (R„ + H.) In 

V„2 y v t stB a 8 t fk T 

x n t ' (R. ♦ HJ +V« + 2R H. + R 

_ 8 l { S l s 


O o o o 

where AR. - -HN E cscE ctnE ■ - AR. E ctnE ■ -HAE. E cscE ■ firat order correction 
to s to to 
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refraction correction vanishes. For the overhead pass assumed here, PC>\ occurs 

o 

at zenith and all the E and R refraction corrections calculated from the formulas 

listed here vanish as they should at E = 90°. 

Table 4 
o 

Limiting Values for E, R, and R Refraction Corrections 
(Using N g = 0.000313, H = 6951.25 Unless Built-In) 


Refraction 

Formula 

Horizon Case (E = 

0°, 

Zenith Case (E = 90°)* 

AE 

(degrees) 

AR 

(meters) 

AR 

(cm/ sec) 

AR 

(meters) 

DC 

00 

98.6 

331.1 

2.74 

Freeman 

00 

00 

00 

2. 18 

NONAME 

1.09 

101.5 

364.7 

2. 57 

GDAP 

0.53 

67. 1 

93.2 

2.25 

NAP-1 

0.63 

86.8 

128. 1. 

2.74 

Laser 

— 

00 

— 

2. 10 

SECOR 

— 

114.4 

— 

2.70 

C-band 

00 

00 

— 

2.22 

Hopfield 

— 

— 

186.5 

— 

Astro 

— 

— 

00 

— 

CRPL 

0.72 

104.2 

— 

— 

REEK 

0.57 

77. 1 

— 

2. 18 


o 

In the Freeman correction formulas for R and R, the restraining function 


does not prevent the correction from becoming infinite at the horizon (E = 0°). How- 
ever, it does serve to improve the nominal correction values for higher elevation 
angles (E >8°), as will be seen later. Dr. Freeman claims validity for his equations 
only for E> 30°. 

GRAPHS OF TROPOSPHERIC REFRACTION CORRECTIONS FOR E, R. AND R 

o 

In Figures 1, 2, 3, and 4. the E, R, and R refraction corrections from 
the different formulas and from the REEK ray traces and CRPL tables are plotted 
against elevation angle. In all cases the tracking station is assumed to be at sea 
ievel. For the radio frequency systems, a refractivity of N g = 0.000313 and a com- 
patible scale height of H = 6951 meters from the CRPL tables are assumed if not 
otherwise specified. The value N g = 0.000313 is the long-term average value at the 

earth’s surface for the United States (Reference 13). For comparison of the optical 

o o 

*For the zenith case (E = 90 \, all formulas for AE and AR yield a zero correction. 
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ELEVATION REFRACTION ERROR ( DEGREES) 



Figure 1. Elevation Error Due to Tropospheric Refraction 
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Figure 2. Range Error Due to Tropospheric Refraction 
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Figure 4. Range Error Due to Tropospheric Refraction, Laser Systems 
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frequency laser formula, an N g = 0.000289 and a compatible scale height of II = 7367 
meters are used in REEK, since they are closer to the average values for the laser 
optical frequency. 

o 

In Figures 5, 6, and 7, the E, R, and R refraction corrections from the 

different formulas are plotted against N for several different elevation angles to show 

the amount of variation which can arise from different values of N . For each value 

o s 

of N in the Freeman formulas forAR anc'AR, a compatible value of H is deter- 
mined from the CRPL tables. A nominal H is built into all the other formulas for 
o 

AR andAR. 

In Figures 8, 9, and 10, the restraining functions which modify the nom- 

o 

inal refraction corrections for E, R, and R are plotted against elevation angle to show 
more clearly the essential differences between the various refraction correction for- 
mulas. 

DEVIATION AMONG REFRACTION CORRECTIONS FOR CONSTANT N 
s_ 

For elevation angles above about 10°, where most satellite tracking is 

done, the refraction error curves in Figures 1, 2, 3, and 4 deviate from each other 

much less than at lower elevation angles. Also, the magnitude of the corrections is 

much less at the higher elevation angles. The deviations among correction formulas 

are shown more clearly in the restraining function plots of Figures 8, 9, and 10. For 

comparison purposes, an effective restraining function for each REEK ray trace is 

obtained by dividing the complete REEK correction at each L by the appropriate 

nominal correction at that E. 

The maximum deviations among the corrections in Figures 1, 2, and 3 
are given in Table 5 for several angles. 

Table 5 

Maximum Deviation Among Refraction Corrections 


E (degrees) 

A E (arc sec) 

AR (meters) 

o 

AR (cm/sec) 

2. 5 

450.4 

33.2 

75.2 

4.0 

162.0 

12. 2 

33. 8 

10.0 

46. 1 

3.5 

2.8 

15.0 

28.5 

2.3 

1.4 

40.0 

9. 1 

1.0 

0.4 


It is evident that some of the spread among the computed refraction 


corrections could be reduced by a different choice of built-in constants in the formulas. 
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Figure 5. Elevation Angle Error versus Surface Value of Refractivity 
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Figure 7. Range Rate Error versus Surface Value of Refractivitv 
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RESTRAINING FUNCTION FOR ELEVATION ANGLE REFRACTION ERROR 
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Figure 8. Ratio of AE^ to AE^q for Elevation Tropospheric Error 
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Figure 9. Ratio of AR t to AR to for Range Tropospheric Error 
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Figure 10. 


Ratio of ARj to AR lo for Range Rate Tropospheric Error 
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VARIATION IN REFRACTION CORRECTIONS WITH CHANGE IN N T s 



In Table 1, all the refraction correction formulas for E are linear func- 
tions of N for a given E, except the NAP-1 formula which uses a constant built-in 

value of N = 0.000350. 
s 

Since the NAP-1 formula is otherwise practically the same as the GDAP 

formula, the NAP-1 formula yields the same correction for E as the GDAP formula 

near N = 0. 000350, but diverges for other values of N . 
s s 

In Tables 2 and 3, the DC, NONAME, GDAP. and C-band refraction 
o 

correction formulas for R and R are ail linear functions of N for constant E. These 

s 

formulas contain nominal built-in values of H = 8750, 8432.336, 7200. and 7600 meters, 
respectively. 


The NAP-1, laser, SECOR, and Astro refraction correction formulas 
o 

for R or R all contain nominal built-in values for the product HN g . These are 

2.7432 meters for NAP-1, 2. 1 meters for laser, 2.7 meters for SECOR, and 2.3 

meters for Astro. Therefore, the NAP-1, laser, SECOR, and Astro refraction 

o 

correction formulas for R and R do not vary with N g . Again, since the NAP-1 and 

GDAP formulas are practically the same, the NAP-1 formula yields the same cor- 

o 2 7432 

reetions for R and R as the GDAP formula near N = = 0.000381. 

s 7200 

o 

The Hopfield refraction correction iormula for R is a linear function of 
N g for a given E, but uses a tropospheric scale height associated with a quadratic 
model of the troposphere rather than the H used with exponential models. 

o 

As shown in Figures 6 and 7, the refraction corrections for R and R 
from the REEK ray-trace and the Freeman formulas are the only nonlinear functions 
of N , since these arc the only functions which introduce a change in H, as indicated 

S q 

in the CRPL tables, whenever N is changed. In general, for E>8 and for normal 

s -6 ° 
variation in N from about 253 to 377 x 10 the Freeman corrections for R and R 
s 

agree best with the REEK corrections. 

As shown in the CRPL tables for some sample tropospheres, changes in 
N arc partialh offset by opposite changes in H, so the product HN tends to remain 

S S o 

constant. The product HN g is contained in all the formulas for R or R refraction 
corrections. However, those formulas which allow N to change, but use a built-in 
value of H, introduce more error than the formulas which use a built-in constant 
value for the product HN . For example, from reference (7) the variation in N 

^ S S 

is about * 60 x 10 about a mean value of 313 x 10 , a variation of about + 20%. 
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From the CRPL model this leads to a variation in HN of less than 8%, as shown in 
Table 6. 


Table 6 

HN Variations from CRPL Tables 
s 


N . < 1# ' 8 > 

%A(N g ) 

H (meters) 

HN g (meters) 

% A(HN g ) 

252.9 

-19.2% 

7920. 85 

2. 00318 

-7 . 9 c /r 

313.0 

0 

6951.25 

2. 17574 

0 

377.2 

20.5% 

5772.81 

2. 17750 

+0. 1% 


Thus, if the center nominal constant value for HN is used in the refraction correction 
o s 

formulas for R and R, normal tropospheric variations may lead to an error of -7. 9% 
for low values of N g , and negligible error for high values of N g . Correcting for N g 
alone without a corresponding change in H leads to overcorrection errors up to 20% 
for high values of N g and up to 11% for low values of N g . 

RELATIONSHIPS BETWEEN REFRACTION CORRECTIONS FOR E, R, AND R 
NOMINAL CORRECTIONS 

Several relationships connecting the nominal refraction correction formu- 
o 

las for E, R, and R are apparent by inspection of these formulas. They are: 

AR 

A E = N ctnE = — rr^ cosE 
to s H 

AR. = HN cscE = HAE. secE 
to s to 

. o o o o 

A R. = -HN E cscE ctnE = - AR, E ctnE = -HAE. E cscE 
to s to to 

d(AR to ) d(AR to ) 

dt dE 


AE ANDAR RELATIONSHIP 

It is of interest to determine whether relationships similar to the above 

also hold for the complete GEOS tropospheric refraction correction formulas for 
o 

E, R, and R in Tables 1, 2, and 3, including those formulas with restraining 
functions. 
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For the Cl DAP and NAP-1 GEOS tropospheric refraction correction 
formulas in Tables 1 and 2, a similar relationship is found to exist: i.e. , 

AE =^^-cosE holds for GDAP and NAP-1 
H 

However, this relationship does not hold for the DC, Freeman, and NONAME formu- 
las. The DC and Freeman formulas for AE are merely the nominal correction 
formulas AE , whereas the DC and Freeman formulas for AR represent a higher- 
order approximation to the true values. It is probably necessary that the level of 
approximation in the AE and AR formulas be the same if the above relationship is 
to hold. The NONAME formulas for AE and AR may also represent slightly different 
levels of approximation, 
o 

AR ANDAR RELATIONSHIP 

It can easily be shown that, for all the GEOS tropospheric refraction cor- 
rection formulas in Tables 2 and 3, the range rate refraction correction formula for 
o 

AR is simply the time derivative of the formula for AR. 

Thus, 

° d (AR) _ ° d (AR) 

AR ‘ lit L_ dE“ 

holds for the DC, Freeman, NONAME, GDAP, and NAP-1 formulas. This is not 
generally stated in the source documents. 

The GEOS range rate measurements are all basically derived by deter- 
mining the number of doppler cycles in a short time interval. The number of doppler 
cycles in this time interval corresponds to the range shift in spacecraft position rela- 
tive to the station during the same time interval, where range shift is measured in 
wavelengths of the transmitter frequency. The measured range shift includes the 
differential effect of the atmosphere on the doppler cycles at the beginning and end 
points in spacecraft position corresponding to the beginning and end points of the 
short-measurement time interval. Therefore, it is physically reasonable that the 
range rate refraction correction should be expressed as a differential in range 
refraction corrections per corresponding time differential or as a derivative, with 
respect to time, of the range refraction correction formula. 

C ONCLUSIONS 

The purpose of this study was to compare the different tropospheric 

o 

refraction correction techniques used with GEOS E, R, and R data to determine 
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whether significant deviations exist between these techniques. If deviations exist, it 
is difficult to know which technique gives the best results, since no perfect model of 
the troposphere is available. However, REEK ray-trace results and the CRPL tab- 
ulated results were included as references for comparison of the results from the 
various formulas, and should indicate which formulas are best. 

The following conclusions summarize the results of this study. 

1. The GEOS tropospheric refraction correction formulas and the 

o 

REEK ray-trace results for E, R, and R are in poor agreement 
for E < 10°. 

2. These formulas and ray-trace results for E >10°, where most 

o 

satellite tracking is done, differ in E, R, and R by insignificant 
amounts in comparison to most operational requirements. How- 
ever, for the GEOS Observation Systems Intercomparison Investiga- 
tion (GOSII), these differences (see page 11) are equivalent in size 
to observed differences in system biases for some of the best 
electronic geodetic tracking systems. Therefore, the systematic 
differences in refraction correction techniques are affecting the 
GOSII results. 

o 

3. The refraction correction formulas for R and R contain the product 
HN g . For best results, whenever N g is changed, the value of H 
should also be changed, using the CRPL model. Otherwise, using 
HN g = constant is preferable to using a measured N g and a con- 
stant H. 

4. For all the GEOS tropospheric refraction correction formulas ex- 
amined in this report, the range rate refraction correction forAR 
is simply the time derivative of the formula for AR. In the source 
documents, the A$ and AR formulas are derived independently. 

5. The Minitrack Interferometer direction cosine tropospheric re- 
fraction correction formulas, which are documented in Reference 2. 
are the proper nominal corrections for angle trackers which 
determine the direction of the incident ray. However, lor inter- 
ferometer systems, the refraction error is much smaller, and it is 
better to make no refraction correction at all to the dual ray inter- 
ferometer data than to use the single ray angle tracker refraction 
correction (see Appendixes A and C). 
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APPENDIX A 

REFRACTION CORRECTION FORMULATIONS 



APPENDIX A 

REFRACTION CORRECTION FORMULATIONS 


Different groups using data from the GEOS geodetic tracking systems 
employ different refraction correction equations. This appendix lists some of the 
different refraction correction equations, and describes how they are being used. 
For comparison of the equations, the original forms and notations are converted 
into a common form and notation. 


A. 1 GSFC ORBIT DIFFERENTIAL CORRECTION PROGRAM (DC) 

The GSFC Orbit Differential Correction Program (DC) was developed 
by the Advanced Orbit Programming Branch at Goddard. It is used for operational 
orbit updating and contains optional refraction correction formulas for interferometer 
direction cosines and for range and range rate. It is documented in reference 2, 
program no. F116. 


A. 1. 1 DC MINITRACK DIRECTION COSINE CORRECTION OR ELEVATION 
ANGLE CORRECTION 


The direction cosine corrections are: 

J Q ~ and m c = equation (2), F116 

where 

Jl , m = measured direction cosines 

/ , m = refraction corrected direction cosines 
c c 

1 + N 

q 1 + A g -N 1 (h,t)/h 

A =N = ( /a - 1) = tropospheric refractivity at the station; a 
s s 

constant for a given station for a given month. 
For example, a typical N g = 0.000313. 

1 + N =1 
o 


A-l 



Keeping the tropospheric part 


m = m (1 + A ) 
c s 


a; = i-; = -Ai 

c s 


A m = m - m = -A m 
c s 


are the changes in observed direction cosines due to refraction in the DC program. 

As shown below, these equations are equivalent to a single refraction 
correction equation for elevation angle. 


By definition 

Z = cosE sinA m = cosE cosA 

where E = apparent elevation angle of incident ray 
A = apparent azimuth angle of incident ray 

Refraction causes small changes, A E, in elevation angle, but in 
azimuth, A A = 0 


therefore, taking the derivatives 
A^ = -AE sinE sinA 


Am = -AE sinE cosA, 


AE = - 


sinE sinA sinE sinA 


A E = - 


sinE cosA sinE cosA 


= A ctnE; 
s 


= A ctnE 
s 


AE = A g ctnE = N g ctnE = AE^ q (radians) (Al) 

This equation is the nominal refraction correction for angle trackers measuring 
the elevation angle of the incident ray and is treated as such for the comparisons 
in this paper. For further discussion of interferometer refraction, see Appendix C. 


o 


A . 1 . 2 


DC' RANGE CORRECTION 


where R 
R 


where 


From equation (12), FI 16 

A R = R- R = C,A csc0. + B N (h, t) csc0 o E -10° 

1 c 1 s 1 s l 2 


From equation (14), FI 16 


AR n = R - R = C, A cscE + B N. (h, t) cscE 
2 c is si 


E> 10 


= measured range 
= refraction corrected range 


For the troposphere alone 


A R = C. A csc0, 
1 Is 1 


E -10 


= cos * (cos0, cosE) 
lm 


0 , = tan 1 T 1 

lm lm 


T t = 0.027759 
lm 


sin0 


tan0. = 0. 027759 * 

lm cos 0 


lm 


lm 


51110 . 

COS ^lm = tanS — ^ 0.999614 from trigonometric tables 

m ^lm 


cos0^ = 0. 999614 cosE 


csc0, = 


1 (1 - cos 2 0 1 )'/ 2, -[l - (0.999614 cosE) 2 ] 


2n'/z 


C^ = 8750 meters 


A -3 



8750 A 


8750 N 


AR = 


1 (l - 0.999G14 2 cos 2 e) I/ 2 (l - 0. 


s 

999228 cos 2 e ) I/2 


(A2a) 


0 


= (HN g cscE) 


8750 


sinE 


L H (l - 0. 


2 2 A I/2 


999228 cos E 


J 


= AR 


~8750 

1 "] 

L h ( 

* 9 \ 1/2 

1 - 0.000772 ctn EV _ 


(meters) E < 10 


R 0 = C. A cscE = 8750 N cscE 
2 1s s 


(A2b) 


= (HN g cscE) = A R q (pp)’ (meters) 

A. 1.3 DC RANGE RATE CORRECTION 

From equation (13) F116, for E < 10° 


o oo r 
ARj = R-R cl « B g 


r(t) • _r*(t) 


dNj (h, t) 

dh - CSC0 2 


2 ° 
-C,A esc 0, ctn0, cos0, EsinE 
Is 1 1 lm 


E > 10 v 


- B s N i ( h > l ) csc 0 2 Ctn0 2 VR 


R \ o 

Tif— (T)) EsinE / 

ms' 7 J 


From equation (15) FI 16, for E>10 


o oo fro —\ 

AR = R-R c2 = S pgK 1 ) ‘ £*<t)J 

-[ C l A s + B s N l< h '‘>] B 


dN x (h, t) 


dh 

ctnEl cscE 


o 


A -4 


where R = measured range rate 

R c = refraction corrected range rate 

The tropospheric part is 


ARj = -C A csc^ ctn0 1 cos0 lm E sinE 


E - 10° 


AR = -8750A (.999614) E sinE esc 0 cos0 

IS XI 


E - 10 ° (A3a) 


-8750N, (.999614) E sinE cosE 

il 3// 2 

1 - (0. 999614 cosE)_ 


- *- HN S E <*" E cscE > PF <- 999614 > 2 « 1. 0. 000772 ctn ^ E l V ] 


o 18743.25 


= I h 

to L H 


(1 + 0.000772 ctn 2 E) 3/ 2 J 


(meters/sec) 



o o 

-C.A E ctnE cscE E>10 

1 s 

o 

-8750N E ctnE cscE E>10 

s 

(-HN g E ctnE cscE)(^j^)= R tQ (meters/sec) 


(A3’o) 


A. 2 FREEMAN RANGE CORRECTION 

The following formulas were formulated by J. J. Freeman Associates, 
Inc. , under Contract NAS5-9782 for the Operations Evaluation Branch (OEBl at 
Goddard and are documented in Reference 5. The puroose of the contract was to 
develop better refraction correction equations for use in the calibration of angle, 
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range, and range rate systems. These formulas are being used by the OEB for 
Goddard Range and Range Rate System (GRARR) quality assurance, GEOS-GRARR 
data validation, and some GEOS data intercomparisons. 


A. 2. 1 FREEMAN ELEVATION ANGLE CORRECTION 


The Freeman elevation angle correction (flat-earth approximation) is 


AE 



(A13) 
ref. 5 


where E = observed elevation angle 

E refraction corrected elevation angle 
c 

H/h is neglected since the scale height H~7,225 meters and target height for 
GEOS, h>10 6 meters pp ~ .007 « 1 


The tropospheric formula is 

AE = N g ctnE =AE to (radians) (A4) 


A. 2. 2 


FREEMAN RANGE CORRECTION 

The Freeman range correction (quasi-flat-earth approximation) is 


AR = R - R 

c 



. 2 

tan a 

R 

s 



(A25) 
ref. 5 


where R 


R 

c 

h 

h' 

R 

s 

a 


N 


elcctromagnetically determined range 
true refraction corrected range (meters) 
satellite height ~ 10'^ km for GEOS 
distance from earth surface to ray path 
distance from earth center to surface 
(ir/2 - E) = apparent zenith angle 

(absolute value of tropospheric or ionospheric refractivity = j ft - l| ) 

L h 2 h 

/|N|dh-^LE /| N | hdh 

o s o _ 
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In the previous equation, if | X | assumes the value 


N N c- h 11 

s 

N surface rcfractivitv 
s 

II scale height ~ 7 km 
- 1000 

h' h, | N | = N s e 7 ~0, 

then the tropospheric formula is 

2 

AH N H(1 - cscE 

s n 

s 

‘ AR o b ■ r 

s 

A. 2. 3 FREEMAN RANGE RATE CORRECTJO ‘ 
From equation (A51), reference 5 

A R = N { V D + V 0 tan a ( 1 - =-11—) ) - 
p l R 6 V R cos / J 

S 


Vg tan a 


/ l«l [>-t (* — ) ] 


where \'g - RE 


3 

As h’ approaches the GEOS satellite height of approximately 10 km, the 
tropospheric rcfractivitv at the target, N jt approaches zero. Therefore the range 
rate correction equation reduces to 


o 

AR = - 


o 

E cosE 
sin 2 E 


| N|dh* 


1 - 3/sin 2 E , h 


h “1 

J |N|h’dh’ 


In the previous equation, if |n| assumes the value 

M-v' hV “ 


( 



then the tropospheric formula is 


o 

A R = 


-HN 

s 


o 

E cosE 
sin“E 



3 

sin E 



o 

= -HN EctnEcscE 
s 


l"l+ — 

l R 


(1-3 esc E) | 


(AG) 


- AR «, [ 1 + nb (1 - 3 

s 


9 

csc“E) 


A. 3 GSFC NONAME PROGRAM 

The NONA ME program was jointly developed by the Mission Trajectory 
Determination Branch (MTDB) at GSFC and the Wolf Research and Development 
Company under NASA contract NAS 5-9756-710. It was designed to produce accurate 
satellite orbits from geodetic quality data. The MTDB and OEB groups at GSFC have 
used this program for both long-arc and short-arc GEOS data intercomparisons. The 
program contains refraction corrections for elevation angle, range, and range rate. 
The current corrections are somewhat more sophisticated than those used in earlier 
versions of this program as reported in reference 24. 


The following formulas were obtained from Reference 19. 
N 


A E = 


2 = v ( »tnF 

(0. 0164 + 0. 93 tanE) s 


.( 0 . 


0164 ctnE + 0.93)1 


(A7) 


= AE 

o 


1 

(0.93 + 0.0164 ctnE) 


+8432.336 N 

AR = -TTTTTT = HN cscE 

(0.026 + sinE) s 


.H 


8432.336 


(1 + 0.026 cscE). 


(A 8) 


A p / 8432.336\ 1 

to \ H ) 1+0.02. 


026 cscE 
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(A9) 


o 

AR 


o 


8432.336 N E cosE 

s 

2 

(0.026 + sinE) 


o 

= HN E ctnE cscE 
s 


~/8432. 336^ 

, 2 *1 
^ sin E 

l H , 

/ 2 
(0.026 + sinE) _ 


8432.336 1 

H ( 1 + 0.026 cscE) 2 _ 

A. 4 GSEC GDAP AND NAP-1 PROGRAMS 

The Geodetic Data Adjustment Program (GDAP) was developed for the 
OEB at GSFC by DBA Systems, Inc. , under contract NAS 5-9860. It is a short-arc 
orbit determination program with error model regression capability and has been 
used on the CDC-3200 computer for GEOS data intercomparisons over data spans of 
about 1/4 orbit. The Network Adjustment Program (NAP-1) was later developed for 
the Mission and Systems Analysis Branch (MSAB) at GSFC by the same company under 
contract NAS 5-10588. NAP-1 extends the GDAP capability to long arcs and to opera- 
tion on the IBM 360 computers at GSFC. Both GDAP and NAP-1 contain refraction 
correction equations for elevation angle, range, and range rate. The functional forms 
of the refraction corrections are the same in both programs, but since different values 
of the constants were chosen, both sets of equations are given here. 

A. 4. 1 ELEVATION ANGLE REFRACTION CORRECTION 

The GDAP, elevation angle refraction correction as given ir. reference 27, 
equation (13. i) , is 




where AZ = Z - Z = -(E - E ) = -Ae 
o c ' c 


Z q = (tt / 2) - E = observed zenith angle 

Z c = (tt / 2) - E c = refraction corrected zenith angle 
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a - 2a = 2N 
1 s 

r = R = earth radius at station 
o s 

a = n -1 = N = refraetivity at station 
o s 

L 1 = 4H/R = = 00. 0045154 

s 63 i 81bb 

H tropospheric scale height (nominal 7200 meters 
assumed in GDAP) 

= radius of earih at observer (0.378,166 meters at sea level) 


AE = 


2N cscE 
s 


sinE +^/sin 2 E + 4H/R ( 


= (N g ctnE) 


■i/l +4H/R, 


2 

esc E) 


(A 10a) 


= Ae 


to\ 
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+ 0 . 


(radians) GDAP 


00045154 esc E 


The NAP-1 elevation angle refraction correction (given in a private 
communication) is 


AE = 


or j cscE 


sinE +/JL + sin“E 


where 


a ] = 0.0007 


9 = 4H/R = 0.004 

** s 


A E 


0.0007 ctnE 


1 + /I + 0.004 esc E 


— = (N ctnE) 

T s 


0.00035 


N 


S 1 +/l + 0.004csc 2 E J 


(A 10b) 


A E 


to 


10035 ( 

*s \ 


0.00035 
N 


0. 004 esc E 


O' 


(radians) NAP-1 
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A. 1.2 


RAW IK UK FRACTION CORRECTION 


The GDAP range retraction correction, from reference 4, appendix C, 
section la, is 


AR. 


2(n - 1) H 
o o 


411 

. , . 2 „ o 

sinL «■ A / sin E + — 

o 


where R { AR - (observed range - corrected range) 

n - 1 + N ground index of refraction 
o s 65 

H q = H = scale height (nominal 7200 meters assumed in GDAP) 

R = R = radius of earth at observer (6,378.166 meters at sea level) 
os ' 

E = elevation angle 


A R = 


2N H 
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. „ , . 2_ 4H 

sinE -i /sin E + tt— 
K 
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= (HN cscE) f Lll£2! 2 

l , . / . (4) 7200 2 


1 V 1 + csc E 

v s 
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= AR 


to 
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, / 28800 2_ 

1 v 1 + — R CSC E 
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In GDAP, R , and \ T are calculated for each pass. For these comparisons, 
^ ^ —6 

a sea level R of 6,378,166 meters and an N of 313 (10 ) are assumed, 

s s 


A R = AR 


to 


m- 


1 + /l + 0.0045154 csc E 


(meters) GDAP (Alla) 
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The NAP-1 range refraction correction is 


AR = 


sinE + Jl 1 sin“E 


where a., 5. 4864 meters 


/ “ = 4H/R -0.004 
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5. 4864 cscE 


= (HN g cscE) 


5.4864 


1 + /I + 0.004 esc E 


HN 


(l +^/l + 0.004 csc 2 E^ 


(A 1 lb) 


= A 


R f 2 -' 74:j2 2 1 

to HN s / 2 

1 +/1 + 0.004 esc EJ 


(meters) NAP-1 


Instead of assuming a nominal H = 7200 meters and calculating N g 

for each pass as in GDAP, NAP-1 assumes a nominal IIN = 2.7432 and a nominal 

s 

I4N = 2. 7432 and a nominal 4H/R = 0.004. 
s s 

A. 4. 3 RANGE RATE REFRACTION CORRECTION 

The GDAP range rate correction, from reference 4, appendix C. 
section 2a, is 


o -AR t E cscE 


A R = 


-2N H 
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E cosE 


^/sin 2 E + 4 H q /R o (sinE +^/sin 2 E + 4H/R g ) A /sin 2 E + 4HR £ 


(A 12a) 
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H(l + / 1 + 0.0045154 csc 2 e\/i + 0.0045154 csc 2 E 


o 
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II / 


^1 +J 1 + 0.0045154 esc 2 e)ViT^ 


0U45154 csc“E 
GDAP 


(— ) 

\sec/ 
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The NAP-1 range rate refraction correction is 


^ r = ~ A R E cosE 
sin 2 E 


-5.4864 cscE 


o 

E cscE 


"■» sin 2 E (l + 0.004 esc 2 e) v 4 - 2 + sin 2 E 


(A12b) 


5.4864 (HN s E ctnE cscE) 


HN 


(l +^1 + 0.004 csc 2 e) v ^T~0. 


004 csc 2 E 


= -Ar 


to 


'/ 2.7432X 

\ HN / 
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( x v ^ 7 0.004 csc 2 E^/l 
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004 esc E^/l + 0. 004 esc E 
GSFC and SAO LASER RANGE REFRACTION CORRECTION 


■1 (—) N. 

\ sec / 


NAP-1 


The GSFC and Smithsonian Astrophysical Observatory (SAO) lasers employ 
pulsed ruby rods operated at an optical frequency with a wavelength of 6943 angstroms. 
At this frequency tropospheric refractivity is primarily a function of pressure and 
temperature with very little dependence on relative humidity. So far. laser data have 
been acquired only at high elevation angles (E >15°), due to transmitter power and 
reflector limitations. Consequently, both GSFC and SAO use a nominal laser range 
refraction correction. From references 17, 18. and 28 this is 

A R = R - R =2.1 cscE = (HN cscE) ( 1777 —) (A13a) 

c s \ hin^ • 


=AR 


to 


(—) 

\ HN s 


(meters) 


and SAO changed in May 1968 to 


AR = 


2. 238 + 0. 0414 PT _1 - 0. 238 h 
sinE + 10 -3 ctnE 


2.238 + 533.5 N 
s 

sinE + 10 3 ctnE 


(A13b^ 


where 


P = station pressure in millibars 

T = station temperature in degrees Kelvin 

N = 77.6 £ x 10 _f> = 0.000289 for nominal P and T 
s I 

h g = station height above mean sea level, assumed to l>e zero. 
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A. 6 


AMS AND ETL SECOR RANGE REFRACTION CORRECTION 


The Army Map Service (AMS) tropospheric range refraction correction 
for SECOR was obtained from reference 20. page 37. 

-ZR) 


AR = 


Kj (1 - e 


sinE + KgCosE 


u 


where 


E = elevation angle 


Kj = 2.7 meters = zenith refraction value 


Kg = 0.0236 = "control constant" 


7 = 


H 7000 


meters 


It= slant range (meters) >10’ meters for GEOS 

e e ~0 for GEOS. This term is needed when tracking 

targets within the troposphere such as aircraft. For any 
satellite tracks it may be ignored. 


n 


A R = — 


2.7 CSCE 2.7 <HN s cscE) 


+ 0. 0236 ctnE 


HN 1 + 0.023 ctnE 
s 


(A 14) 


AR to HN, (l + 0.023 ctnE) ( meters ) 

The Army Engineer Topographic Laboratory (ETL) Tropospheric range 
refraction correction for SECOR is also given in reference 20. page 49. This is the 
same as that given above for GDAP in equation (Alla). 
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APL TRANET RANGE RATE REFRACTION CORRECTION 


The Applied Physics Laboratory (APL) tropospheric range rate refraction 
correction for TRANET was obtained from reference 8, equations 8 and 9. 

A R = -N g R s E [f(E)J = - (HN, E ctnE cscE) -j|- sin 2 E [f (E)] (A15) 

o R or i 

* - AR to^f s '" e M 


H 


o 
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( 



R= 23,000 meters = best fitting equivalent thickness of troposphere 
above station 

R g = 6,378,163 meters in this comparison 


This formula was derived from a quadratic rather than an exponential fit 
to the tropospheric refractivity profile. Consequently it looks quite different than 
some of the other range rate refraction formulas. 

v 

More recently, in reference 25, the same author has further developed 
and refined this technique for range rate refraction correction by using two separate 
quartic fits to the wet and dry components of the troposphere. 

A. 8 NWL TRANET RANGE RATE REFRACTION CORRECTION 

The Naval Weapons Laboratory (NWL) tropospheric range rate refraction 
correction for TRANET is given in equation (A6) of reference 21 and in reference 26 
describing the NWL ASTRO computer program. 

Af - 4 s £[•><" -!>]=-! A R 

cos Z 

where Af = refraction correction to Doppler frequency 

. o 

AR = range rate refraction correction 
Z = zenith angle = ( tt/ 2 - E) 
f = satellite transmitter frequency 

Q 
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So 


A. 9 


A. 9. 


A. 9. 


n = mean tropospheric refractive index between station and satellite 
h = height of satellite above the station 

(n - 1) = 2.3 meters for a station at mean sea level 


u 


= -^f E E [h <n - 1>] = -2.3 E ctnE CSCE = (A16, 

WALLOPS STATION C-BAND PROGRAM 

The Wallops Station C-band program is documented in reference 29. 

1 WALLOPS C-BAND ELEVATION CORRECTION 

AE = N ctnE = E. (A17) 

s to 

2 WALLOPS C-BAND RANGE CORRECTION 

A R = 7600 N cscE (A18) 

s 



f 1 
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APPENDIX B 

ELEVATION RATE DETERMINATION 
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APPENDIX B 

ELEVATION RATE DETERMINATION 


The required expression for elevation rate is derived from the geometry 
shown in Figure B-l. 


/ 


V 


TRACKING STATION 



Figure B-l. Plane cont lining satellite, trackingstation, and center of earth. 
(A circular overhead orbit around a spherical stationary earth is assumed.) 


. sin cut = RcosE/R,„; eosE = R^sintut/R (Bl) 

cosut = (R s + RsinE)/R T ; sinE = (R T cos cut - R^)/R (B2) 

R = (R” + R™ - 2R_R cas cut) 1/2 (B3) 

s 1 is 

substitute (B3) into (Bl). 

C 

B-l 


cosE = R~ sin cut (R“ + R 2 - 2R_,R cos cut) ^ /2 
I s l I s 

then by differentation, 

d F - 1 - q 9 

-sinE -jjp * u/R.j,H cos cut - cuR^R * (R^,R g sin cut) (B4) 

substitute (B2) into (B4) 

-1 HF -12 2 2 

-R (R T coscut - R g ) = wR T R~ (R coscut - R T R g sin cut)/R 

jy o 9 9 

-r— - = - cuR~(R„,coscut + R“coscut -2R...R cos cut - 
ut 11 s is 

R R + R R COS 2 cut)/R“(R„coscut - R ) 

ISIS 1 s 

= -<uR t (R t - R g Cos cut) (R T coscut - R g )/R 2 (R T cos<ut - R g 

dF 2 2 

= - w R t (R t - R g COs cut) / (R t + R g - 2R T R g cos cut) (B5) 

R t = distance from the center of earth to the satellite (7, 711,499 meters) 

R = distance from the center of earth to the station (6,378, 166 meters) 
s 

R = distance from the station to the satellite 
cu = angular velocity of the satellite = 2 ir/T 

T = period of the satellite = T (R T /R^)*^“ page 96 of reference 22 

= 84.347 minutes, the period of a circular orbit at the earth's surface 
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APPENDIX C 


MINITRACK INTERFEROMETER TROPOSPHERIC 
R E FR A CT 10 X COR R E CT IO N 


The Minitrack interferometer measures the phase difference between 
the signals received by two ground antennas spaced a known distance apart (see 
Figure C-I). If the propagation time of a wavefront from the satellite to antenna Aj 
is tj. and to antenna A. ; is t.,, then the interferometer measures the time difference, 
At - (t 9 -t >. between these signals. 
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Figure C-l. Interferometer Geometry (No Atmosphere) 

The effect of the troposphere is to bend the two ray paths and retard the 
signals along these ray paths from the satellite to the two antennas. Thus, the propa- 
gation time to antenna A^ increases from t^ to t^ +At , and the propagation time to 
antenna A 9 increases from t. ; to t. } -At.,, causing the measured time difference to 
increase from (t.,-^) to (t., +At. ) i - (t^ + 4 - 1 ^), a net increase ofS(At)=(At 9 -At^. 
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In other words, the difference in retardation along each ray path due to the atmos- 
phere eauses a slight increase of 8(At) in the measured At. 

If the two rays were parallel, as from a distant radio star, and it the 
atmosphere consisted of horizontally homogeneous plane parallel layers, the re- 
tardation along the two rays would be identical, so the At. ; =At^, and 8 (At) = 0 . In 
this case, there is no change in the measured At due to the introduction of an atmos- 
phere; and, since the true angle to the star is also unchanged, there is no need to 
make a refraction correction to the derived interferometer angle. It is primarily 
because the two ray paths to a near-earth satellite are not quite parallel and because 
the atmosphere is not horizontally but spherically stratified that a differential 8(At) 
is introduced. 


For the Minitra ck interferometers, the baselines are horizonta 1 at the 
midpoint and the longest oasthne, B, is less than 125 meters: whereas the range. R. 
to a near-earth satellite is at least 150 kilometers. Thus B « R. Under this condition, 
and assuming that no atmosphere is present, the interferometer direction cosine can 
be expressed as follows: 


m 

o 


cos/3 


cAt 

B 


AR 

B 


where: ni Q = direction cosine for no atmosphere 

/3 = interferometer angle between direction line to satellite and baseline at 
baseline center. 


c = velocity of light in vacuum 


B “• length of baseline in meters 

At = t 2 - 1 ^ = difference in propagation time from satellite to antennas A i; and 
A 1 as seen by interferometer. 

AR = R., -R^ = difference in propagation distance (meters) from satellite to 
antennas A 9 and A^ (derived from cAt) 


Also, from Figure C-l, the change in /3, A/3, is 
A/3 = /3 x -/ 3 2 ^ for B «R. 
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When the troposphere is introduceo, the measured propagation time 
difference, At, is increased by the differential retardation 8 (At) = At 0 -At^ so that 
with no correction for refraction, the measurement is interpreted as: 


m - cos /i 


[ At + 8 (At)] c = AR + 8 (AR ) = 

B B m o 


+ Am 


The change in the derived direction cosine due to the troposphere is 


therefore: 


a — m ~ 

A m - m - m - — — — 

o B 


(Cl) 


where: 8 (AR) = differential change due to the troposphere in the derived range 
difference along the two ray paths. 


From Table 2 on page 4, the nominal derived increase in a single path 
range observation due to the tropospheric retardation is given as 


A R. = HN cscE 
to s 


where: H = tropospheric scale height (~7 x 10 meters) 


,-4. 


N g = refractivity at station (~3 x 10 ) 


E = elevation angle of satellite 


Thus, for the interferometer, the effect of the troposphere is to increase 
the derived range difference between the two ray paths by: 


8 (AR) = <AR tQ ) 2 - <AR tQ ) 1 = HN g (esc E 2 - cscE^ 
~ d(^ R to) . 


d E 


AE = (-HN cscE ctnE)AE 
s 


where: 


(ARto >2 = predicted increase in range observation along ray path to antenna 


Ag due to tropospheric retardation. 


(AR tQ )^ = same for the ray path to antenna A^ 


AR. 


d 

■ , ■ - = -HN cscE ctnE = rate of change of predicted tropospheric 
dr s 

increase in derived range with a change in elevation angle, E. 


AE = Eg - Ef = small change in elevation angle between the ray paths to 


the two antennas due to the geometrical separation of the antennas. 
The small correction toAE, due to the different changes in K 
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and E (? caused by the troposphere, is ignored here. For E> 10°, 
the error in A E due to this effect is 


8(AE) 


d(AEl A _ a< N s CtnK ) 

aE AL at 


AE = -N csc“E AE, which is S . 01 AE . 
s 


The difference in elevation angle, AE, for the two ray paths to the two 
antennas is related to the difference in interferometer angle, A /!, for the two ray 
paths to the two antennas as shown in Figure C-2. 
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Figure C-2. Relation Between AE and A/3 

At the satellite, the angle A/3 is defined as the angle subtended by the 
baseline. Therefore, the projection of the baseline on the great circle trace c i the 
interferometer angle (i produces the short arc A/3 at the satellite. TheAE corre- 
sponding to this A/J is also shown in Figure C-2. The arcs^ and — form two 
sides of a small right spherical triangle with included angle a. Since the spherical 
triangle is small, it may be replaced by its projection on the tangent plane as shown 
in T-'jgurc C-2. 

Then AE - A/3 cos a - — cos a. (C2) 
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To determine cosc, return to the large right spherical triangle with 
sides K, p. and A. From spherical trigonometry, 


cos a = tanE ctn/3 
Now returning to equation (Cl) 

Am = - 8( £ R) = 4 (-HN cscE ctnE)AE 
B B s 

Using equation (C2) 

Am = 4 (-HN cscE ctnE) — 5 cos a 

B S n 

Using equation (C3) 


Am = 4 (-HN cscE ctnE) (tanE ctn/3) 

hJ S K 


,-HN cscE. 
= 1 s > 


ar. 


ar. 


n to to 

cos/3 = s— cos/3 = 5 — m 

K K 


or 


A AR„ 

A m to 


m R 

For the flat -earth approximation (see Figure C-2) 


Then 


R = H t cscE 
Am _ AR to 


m 


R 


HN cscE „ 

= - H N 

H t cscE H s 


Or since m = cos/3 

Am = A(cos/3) = -sinjSA/3 
a . Am Am . _ H XT . , 

= - 5553 ' - — c,n(3 = h/ V 1 "'* 


For nominal values of scale height, H = 7 x 10 meters, and of GEOS 

3 

satellite height, H = 10^ meters, = - X - -5* ~ 0.007. 

1 H t 10 6 


This result (C7) differs considerably from that given in reference 2, 


where: m = — = m (1 + A ) 
o q ' s 


or 


A m = m - m = -A m 
o s 

A21=. a = - N 

m s s 

A 3 = ctn/3 = N ctn/3 
m s K 


(C3) 

(C4) 


(C5) 


(C6) 

(C7) 


(C8) 

(C9) 
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The Minitrack interferometer correction for the troposphere derived in 
(C7) says that the correction from reference 2 should be reduced by a factor of 
about 0.007 for the GEOS satellite. 


On the other hand, (C7) is equivalent to that derived for the Minitrack 
interferometer aircraft calibrations in reference 30, where 

m = mu = m (1 + N) 
o ^ 


m ' m o . Am 
m m 


= - N 


or A p = - ctn,3 = N ctn,? 

Using an exponential model for the troposphere, it can be shown that 
HN 


(CIO) 

(Cll) 


N = 


H. 


Assume N(h) - N e 
s 


-h/H 


where it is implicit that h is height above the tracking site, not above sea level, 
and N is the refractivity at the site. 




N(h) dh 


Then N = ° 




N s f H t -h/H HN s r 

‘hJ e dhI T 1 


dh 


t o 


l-e' H * /H ] 


-H t /H -1000/7 
where e = e — 


N 



and (C7) and (Cll) are equivalent for satellite heights. 

To show that (C7) is equivalent to (Cll) for aircraft heights, it is better 
to rederive (C7) in terms of height above sea level instead of height above the site. 
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Then, instead ofAR. = HN cscE, 
’ to s 


AR 


to 


= J 1 N(h) dr =J 1 N q e " h/H dh cscE = 
R H 


-H /H -H./H 


-u / ii -n / 11 “i r *i 

H cscE Ne s - N e 1 = H cscE N - N . 

. o o J L s tj 


where: N y = refractivity at sea level 

N g = refractivity at site 
N t = refractivity at aircraft 


Then 


Am . AR to _ H(N s - N t> cscE . H(N s - N t> 


m 


R 


(H. - H ) cscE H - H 

' t s t s 


(C12) 


* H(N - N ) 

- ~ h. s ~ ' h ^ 

t s 


(C13) 


For this relation to be equivalent to (Cl 1) for aircraft heights, it must 


be shown that N - 


H(N g - N t ) 

H - H 
t s 


J 1 N(h) dh J 1 N Q e “ h/H dh 

H H 


N = 


/ H ‘ 


dh 


H - H 
t s 


-H /H -H /H 

H(N e s - N e 1 ) 

o o 

H - H 
t s 


H(N s - N t ) 
H t" H s 


so (C7) and (Cl 1) are also equivalent for aircraft heights. 


The advantage of this new derivation of (C7) or (C 11) is that more precise 
expressions for the single path range (or range rate) refraction correction than those 
used here are available and can be used for a more precise interferometer refraction 
correction at low elevation angles, subject to the conditions that B« R and the inter- 
ferometer antennas are at the same height. 
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As an example, a more precise expression of 8 (AR) is available from 
the range rate refraction correction equation for GDAP from Table 3 on page 5: 


S(A RI = ^ 


_ d(AR) 


dt 


At = 


(AR^At 


7200 2_ 

H / T 

+ 0. 0045154 CSC E + ( 1 + 0. 0045154 esc E) 


AR (o At f(E) = (-HN g cscE ctnE) AE f(E) 




Returning to (Cl) 

Am = 4 8 (AR) = 4 (-HN escE ctnE)AE f(E) (C14) 

n 13 S 


(C15) 


H(N - N ) 

„ ~ f(E) ctn/J = 
t s 

N f(E) ctn ii (C16) 

A similar derivation to that used to obtain (C15) and (C16) for the 
troposphere should apply to the Minitrack refraction correction for the ionosphere. 
This will be investigated at a later time. 


Except for the factor f(E), (C14) is equivalent to (C4) 
It follows that 
AR. 


and 


f ( E, 

m R 

A/3 = ctn/3 = f (E) ctnp = 


AR 


to 
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